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Abstract 

The possibility of using the optimized 5 expansion for studying medium 
effects on hadronic properties in quark or nuclear matter is investigated. The 
5 expansion is employed to study density effects with two commonly used 
models in hadron and nuclear physics, the Nambu-Jona-Lasinio model for the 
dynamical chiral symmetry breaking and the Walecka model for the equation 
of state of nuclear matter. The results obtained with the 5 expansion are 
compared to those obtained with the traditional Hartree-Fock approximation. 
Perspectives for using the 5 expansion in other field theoretic models in hadron 
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and nuclear physics are discussed. 
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I. INTRODUCTION 



The study of possible modifications of hadron properties in the nuclear medium is one 
of the central problems of contemporary nuclear physics. In principle, these and related 
phenomena in nuclear physics are governed by the fundamental theory of the strong interac- 
tions, quantum chromodynamics (QCD). However, although QCD has been very successful 
in explaining a large class of hadronic processes at high energy and large momentum transfer, 
typical nuclear phenomena at lower energies cannot be derived from QCD with the theoreti- 
cal tools presently available. The difficulty of using QCD for phenomena at the nuclear scale 
is related to the nonperturbative nature of these. Due to the asymptotic freedom property 
of QCD, high energy processes are calculable by perturbative techniques in the quark-gluon 
coupling constant. On the other hand, since there are no reliable systematic approxima- 
tion schemes in field theory for performing nonperturbative calculations, the construction 
of models is an important aspect of low energy QCD. While there is considerable optimism 
that eventually one will be able to solve QCD numerically on the lattice using supercom- 
puters, the development of analytical approximation methods are in urgent need to make 
contact with the wealth of data on nonperturbative phenomena, presently available, or that 
will be available when the new experimental facilities under construction start operating. 
The 5 expansion []!]] is an example of a method recently developed with the aim of studying 
nonperturbative phenomena in field theory. 

The idea of the 5 expansion is to perturb the original theory by the introduction of an 
artificial expansion parameter S, absent in the original theory. The parameter 5 is introduced 
in such a way that it interpolates between the theory one wants to solve and another theory 
that one knows how to solve. The 5 expansion can be formulated in two different forms, the 
logarithmic S expansion [0 and the linear S expansion 0- |3| . In this paper we consider the 
linear form. Specifically, let S be the action of the theory one wants to solve, and So the 
action of the soluble theory. Then, the interpolating action S(S) is defined as 

S(5) = (1 - 5)S + SS, (1) 

so that S(0) = S and S'(l) = S. The next step involves the evaluation of desired physi- 
cal quantities as a perturbation series in powers of 5, which is then set equal to 1 at the 

3 



end. A crucial aspect of the method is the recognition that So involves arbitrary unknown 
(dimensionful and/or dimensionless) parameters. If one were able to solve the new theory 
to all orders in 5, the unknown parameters would not play any role, since no physical ob- 
servable would depend on them. However, since we will be able to solve the interpolated 
theory only to a finite order in 5, there will remain a residual dependence of the results on 
the parameters of Sq. These arbitrary parameters must therefore be determined according 
to some criterion and in fact there are many ways in which this can be done ||- |J. One 
physically appealing way to fix the unknown parameters, which is the one adopted here, is 
the principle of minimal sensitivity (PMS) introduced in Ref. ||. This principle amounts 
to the requirement that physical quantities should be at least locally independent of the 
parameters. In the original applications of the method, the unknown parameters were set 
to be equal to unity. The S expansion, together with the criterion of the PMS of physical 
observables, is known as the optimized S expansion. The convergence of the optimized S 
expansion has been proved in Ref. |7]]. 

The different forms of the S expansion have been successfully applied to many different 
problems in quantum mechanics 0, particle theory |8||J, statistical physics |L0| and lattice 
field theory PjrT|. Motivated by these successes, in this paper we investigate the possibility 
of employing the linear 5 expansion to study medium effects in hadron and nuclear physics 
using typical field theoretic models. In the next section, we consider the use of the 5 
expansion in the study of density effects on the chiral symmetry breaking in the Nambu- 
Jona-Lasinio (NJL) model. The PMS criterion for the typical chiral quantities is applied 
following previous experience with the method. We also investigate an alternative way 
of fixing the parameters by applying the PMS to the energy density of the system. In 



Section |T| we consider the Walecka model and study the effective nucleon mass in nuclear 
matter. Both problems can be treated by obtaining the propagators of the fields involved. 
Traditionally, the propagators are obtained in the Hartree-Fock (HF) approximation. In 
the case of the NJL model this approximation amounts to neglecting the corrections to the 
four-point vertex. In the Walecka model, the HF approximation consists in obtaining self- 
consistently the nucleon propagator with bare meson propagators, and neglecting corrections 
to the meson-nucleon vertices. An additional approximation in the Walecka model is the 
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neglection of vacuum effects in the nucleon propagator. With the purpose of comparing the 5 
expansion method with the traditional HF approximation, we also neglect vertex corrections 
in both models, and neglect the vacuum in the Walecka model. Conclusions and discussions 



of the perspectives for future calculations appear in Section IV. 



II. NAMBU JONA-LASINIO MODEL 



In the limit of zero current quark masses, the two-flavor Lagrangian density of the 
Nambu-Jona-Lasinio (NJL) model ]nj is given by 



£ N jl = 1^0)1 + G (<?<?) 2 - (l^rq) 2 



(2) 



where the quark field operators q = q(x) represent the doublet of u and d quarks. Since the 
model is non-renormalizable, one has a cutoff A as an extra parameter of the model, besides 
G. 

According to Eq. (JJ), one needs to introduce a Lagrangian density C such that 



£njl(<5) = (1- 6)C + 5£ 



NJL 



(3) 



where the equations of motion derived from C can be solved as exactly as possible. Since 
we are looking for solutions which break chiral symmetry, the natural choice for Co is 



Co = q(i0- n)q 



(4) 



where /i is an arbitrary mass parameter introduced for dimensional reasons. Therefore, the 
interpolated NJL Lagrangian density can be written as 



£njl(*) = (1 - 5) [q(i0 - fj,)q] + 6 [q{i$)q + G [{qqf ~ [qisTqf] } 
= q{i$ -fi)q + 6 {G [{qqf - (q^rq) 2 ] + fiqq} . 



(5) 



The evaluation of physical quantities is performed using perturbation theory in the pa- 
rameter 5. The physical quantities of interest, whose values characterize the chiral symmetry 
breaking, are the constituent quark mass M q , the quark condensate < qq > and the pion 
decay constant f n . The quark condensate for a given flavor is given by 
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< QQ >= J (^Tr [S(p)\ , (6) 
where the trace is over spinor and color indices. The pion decay constant can be evaluated 



using the Pagels-Stokar JTJj formula: 

r-b 



iUq„S ab = g nqq J TT^TjTr 



l^T a S{q)^—S{q + p) 



(?) 



d 4 p 
(2?r 

where the trace now is over spinor, flavor and color indices. The quark-pion coupling is 
obtained from the Goldberger-Treiman relation. 

In order to calculate these quantities one needs the quark propagator S(p), which can 
be obtained using Dyson's equation. Expressed in terms of the self-energy S(p) the quark 
propagator reads 

5- 1 (p) = 5 - 1 (p)-E( P ), (8) 
where Sq (p) is the inverse of the quark propagator corresponding to Co- 

So 1 (p) = ff- f i. (9) 

S(p) is calculated as a power series in 5, and S(p) is then obtained by inverting Eq. 
Since the self-energy is calculated perturbatively, this is still a perturbative scheme. It is the 
application of the PMS to physical quantities which furnishes the nonperturbative character 
to the 5 expansion. 

In zeroth order in S, one has 

sW(p) = 0. (10) 

Therefore, one can invert Eq. (|8|) by using the well-known modified Feynman ie prescription 
for particles in a Fermi sea |Tj| . The in- medium quark propagator is then given by: 

- + *i£§f V - *G») (Pr - |p|) , (11) 

where Eq{j>) = (p 2 + /i 2 ) 2 , and Pp is the Fermi momentum, which relates to the density p 
via P F = {n 2 p/2) l /\ 

At this zeroth order in 5, no dynamical content from the model has been used. The 
dynamics of the model starts to show up at order 5. At first order in 5, the self-energy E(p) 
is given by: 



E«(p) = -Sfi 

+ 2i8G J {Tr [5 (ff)l - 5 (g) - 75 r a Tr [r a S (q) l5 ] + i 5 r a S (q)r a ^} , (12) 



where a sum over the isospin index a is implied. 

Substituting Eq. ([11]) into this equation, we obtain for the following expression: 



where 



M x = 5^fi (N c N f + I 



and 



S (1) (j9) = -^ + M 1 -7oS 



A(A 2 + /i 2 ) I -P F (Pl + /U 2 ) 5 -/i 2 ln 



(13) 



2^1 



A + (A 2 + /i : 
P F + (P 2 +/^ 



(14) 



^o = -4^/^(P,-|q|) 



(27T)* 



(15) 



Since the effect of So is just to shift the chemical potential [[L4| one may write the constituent 
quark mass to 0(5) as 



M q = n-dfi + Mt. 



(16) 



Substituting Eq. ( |TTD into Eq. @ and Eq. (|7|), one gets for the order parameter per 
flavor and for the pion decay constant the following lowest order expressions: 



< qq > 



2N c fi 



a(a 2 + - p f (p 2 f + //r - /i 2 in 



A + (A 2 + /x 2 )^ 
P F + (P F + ^ 



(17) 



and 



f 

■ J 7T 



2 _ N C N } ^ 



(2k) 



In 



2\2 



A+(A 2 + /i 2 ) 
P F + (P F + ^ 



(18) 



where the lowest order Goldberger-Treiman relation, g nqq = fi/ f w , has been used. 

The next step in the process is to fix \i. In a previous work |D| this arbitrary parameter was 
determined by requiring /„., which is a fundamental quantity in the study of chiral symmetry 
breaking, to be stationary with respect to variations in /i. This is also a convenient choice 
since, apart from having a well known empirical value, f n is the only one of the studied 



quantities in Ref. || which has a well defined stationary point for finite values of /i. To fix 
the noncovariant cutoff A one uses the empirical value /„. at zero density. If one applies the 
criterion of stationarity to f % in the vacuum, one obtains (for N c = 3 and Nf = 2) 



With the input /„. = 93 MeV, one finds A = 571 MeV, which implies jj = 553 MeV and 
< q q >= -(250 MeV) 3 . Using GA 2 = 2.89 as in Ref. § and setting S = 1 in Eq. pD 
one finds the constituent quark mass to 0(5) to be M q = 574 MeV. Figure 1 shows \i as 
a function of the density and has been obtained by applying the PMS to f n for different 
values of Pp. The results obtained for M q and — < qq > 1 / 3 , for different values of Pp are 
shown in Figure 2 (solid and dashed lines respectively). In Figure 3 the solid line shows the 
Pp dependence of f^. Contrary to what happens at finite temperature || we find that our 
result are insensitive to whether n(Pp) (as in Figure 1) or /x(0) = 533 MeV is used. 

A natural question which arises at this point is the uniqueness of the of the value of /i. 
If one were to use other physical quantities to fix /i, as for example the masses of the vector 
mesons, it is very likely that one would obtain a different value for fi. In such a case, one 
would have a different quark propagator for each observable. Of course this would not be 
catastrophic if the spread of the values of [i determined with different observables is not 
too large. In order to avoid such potential uncertainties, we propose to fix \x by demanding 
that the energy density of the system be stationary with respect to variations of fi. The 
energy density can always be written in terms of the propagators of the theory, and then 
it is natural to demand stationarity of the energy with respect to the unknown parameters 
of the propagators. Then, all physical observables are determined from the same quark 
propagator. 

From the Lagrangian density, Eq. @, we have that the energy- momentum tensor is given 



Note that we have not used the equation of motion for the quark field operator. The energy 
density is the volume integral of the expectation value of T 00 in the many-quark state. The 



H = 0.97 x A . 



(19) 



by: 



t njl = iqr 



»d v q - £njl 




(20) 



8 



expectation value of the field operators can be evaluated using the usual Wick contraction 
technique. This leads to 



£ N jl = y I <r\r < r m > 



—I 



+ 



A Or 



» 4 

d 4 q d A k 
(2tt) 4 (2tt) 4 
d 4 q d k 



g°Tr ^°S(q) 



+ i 



d 4 q 



Tr yS(q)] - G 



d 4 q 



Tr [S(q)} 



Tr [S(q)S(k)} + 



d 4 q 



(2vr) 4 

Tr b 5 T a S(q) l5 r a S(k)}) 



Ti[Tj 5 S(q)} 



(2tt) 4 (2tt) 4 } 
Substituting Eq. (O) into the expression above, we obtain: 



(21) 



(0) 
NJL 



-2N c N f / A f Q F q . - 2GN c Nf{2N c N f + 1) 



d 3 q \i 
p F (2tt) 3 E (q) 



lp F (2n)3E (q) 

The requirement that £ be stationary with respect to variations in \x leads to 

1\ f A d 3 q fi 



(22) 



H = AG[ N c Nf 



(23) 



2) Jp f (2n)3E (q) ' 

This is is the familiar Hartree-Fock gap equation of the model, where \x has the interpretation 
of the dynamically generated mass. 

Next, we consider the first order self-energy. By inverting Dyson's equation, Eq. one 
obtains the quark propagator: 



S (1, (p) 



A + Ml + -*&>))» (ft- -IPI) 



p{ - Mf + ie 



(24) 



where 



rf = (rf,P) = (P° + So,p) 



Ei{p)= P 2 + (M X 



(25) 
(26) 



The superscript (1) in indicates that the propagator has been obtained with a self-energy 
calculated to first order in 8. Note that we are not expanding the propagator in powers of 
5. The process of obtaining the propagator by inverting Dyson's equation with a self-energy 
calculated in perturbation theory is known as the chain approximation. 
Substituting Eq. (|4]) into Eq. (|2l|) , we obtain: 
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f & = -2N C N S £ (05^ - KOT^W/ + 1) 
where Mi is given by Eq. (THI) and Ei(q) is defined in Eq. 



A d 3 q Mi 
p F (2n) 3 E 1 (q) 



(27) 



Application of the PMS to £njl> 

dfx dMi dfj, 

leads to 



(28) 



Again, we have obtained the familiar Hartree-Fock gap equation for the dynamically gener- 
ated mass. 

If one proceeds to higher orders in 5, in the scheme of neglecting vertex corrections, the 
higher order quark propagator will always be of the form of Eq. fl24|), with Mi replaced by 
another constant, say M, which is a function of \x. However, because of the PMS condition 
on M at each order will always be given by the same value. This value is the one that 
satisfies the usual gap equation: 

(30) 

where 

£(g) = (q 2 + M 2 ) f . (31) 

Therefore, the PMS condition on the energy density is equivalent to the usual Hartree- 
Fock solution for the dynamically generated mass, in the approximation of neglecting vertex 
corrections. 

This result should be compared to the one presented in Ref. 0] where, in the context 
of the effective potential, it was found that the 5 expansion and the 1/N expansion are 
identical in the large N limit. 

In Figure 2 we compare the results obtained for the quark mass and the quark condensate, 
when the two above described ways of applying the PMS are used. We call PMS1 the 
results obtained by imposing f„ to be stationary with respect to fi, and PMS2 the results 
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obtained when the PMS is imposed to the energy density. The solid and dashed lines give, 
respectively, M q and — < qq > x / 3 obtained with PMS1, and the dotted and dot-dashed lines 
give, respectively, M q and — < qq > x / 3 obtained with PMS2. In both cases we used the 
same set of parameters: A = 571 MeV and GA 2 = 2.89. 

In Figure 3 we show the results obtained for f n with PMS1 (solid line), PMS2 (dotted 
line), both with the parameters given above, and with PMS2 (dashed line) with a new set 
of parameters: A = 653 MeV and GA 2 = 1.98. The last set of parameters was fixed by 
renormalizing /„- and < qq > at Pp = to their experimental values and requiring M q to 
be roughly one-third of the nucleon mass. The curves for the quark mass and condensate 
obtained by PMS2 with the renormalized parameters are not shown because their behavior 
is analogous to /„•: they go to zero at Pp = 1.6 fm _1 . At Pp = we have M q = 314 MeV. 
These are the usual HF results and the renormalized parameters are the same as used in 



Ref. |B] for two flavors and a three momentum cutoff. 

From these figures we see that the results change appreciably when different criteria are 
used. The main difference is related with the density dependence of the quantities: while 
with PMS1 the quantities smoothly approach zero at some critical density, they go to zero 
through a first order phase transition with PMS2 (or Hartree-Fock) 

III. WALECKA MODEL 



In this section we consider the Walecka model |[5] for nuclear matter. The Lagrangian 
density of the model is given by 

Av = i> - g u V) — (M — gM V> + ~(<^W - m 2 ^ 2 ) - -F^ v + \mlV^ , 

(32) 

where if) represents the nucleon field operators, and are respectively the field operators 
of the scalar and vector mesons, and F^ u = d^V u — d u V^. 

The energy-momentum tensor density corresponding to this Lagrangian density is given 



by: 

T% = nPYd u ip + d"(j)d u (j) + d'VxF^ - g^Cw . (33) 
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Note that we have not used the nucleoli equation of motion. Next, we eliminate the meson 
field operators in favor of the nucleon field operators. The Euler-Lagrange equations yield 
the meson field equations: 



+ ml) V v = g^YiP 



(34) 
(35) 



In obtaining the second equation above we have used baryon current conservation, which 
implies that <9 M V M = 0. These equations can formally be integrated as: 



(j){x) 



-fjo 



d 4 yA a (x - y)ip(y)4>(y) , 
<PyA w (x - y)^(y)j^(y) , 



where Aj(x), i = a, u, is given by: 



AJx) 



d\ 



-iqx 



(36) 
(37) 

(38) 



(27r) 4 q 2 - mf + ie 

In Eq. (|37|) above, because of baryon current conservation, we have neglected the term 
proportional to p^p" /m^ in the vector meson propagator. 

Using the expressions above for and in Eq. (^), taking the expectation value of the 
resulting expression in the many-nucleon state, and evaluating this with the help of Wick's 
contraction technique, we obtain: 

d 4 q 



with 



2 777.2 
1 



(2tt) 4 
d A q 



{Tr [ftf - g 



M)]%)}+(n,+(Ta, 



(39) 



■TrS{q) 



(2vr) 4 

(q 2 - ml) AM 2 ) ~ 1 



//"" - 9l I tSttSi^ [S(Q + k)S(k)] AM 2 ) 



(2tt) 4 (2tt) 4 

f-?Y\(? 2 )}, 



(40) 



and 



1 q 2 

/rpfJ,V\ Huj 

2 m z 



d A q 



+ 9i 



(2tt) 4 

d^q d 4 k 
(2vr) 4 (2tt) 



Tr7 M 5(g) 
Tr f 



d\ 



Tr 7 ^(g) 



(2vr) 4 

lx S(q + k)j x S(k)\ A^q 2 ] 



fJ,U 



X 



■{q 2 -ml)A^)-l 



(41) 



12 



In the same way as in the NJL model, the nucleon propagator is obtained by inverting 
Dyson's equation 



S-\p) = S 1 (p)-Z(p), 
where S is the propagator corresponding to £ i n Eq. (|55| ) below, 

S \p) = p / -M , 



(42) 



(43) 



with the self-energy S(p) calculated as a perturbation expansion in powers of 5. 

For infinite nuclear matter, because of the translational, rotational, parity and time 
reversal invariances, S(p) can be generally written in terms of the unit matrix and the Dirac 
7^ matrices as follows [ p~6f| : 



S(p) = S s (p)- 7 ^(p) 

= £ V, |p|) - 7 °£ V, IPl) + TP^V, |p|) . 



(44) 



Defining the following auxiliary quantities flq| : 



m*( p ) = m + j: s (p) , 

p* = p[l + E»] , 
E*{p) = [p* 2 + M* 2 {p)f , 

p*» = p » + 57(p) = [p° + E°(p), p* 

we can invert Eq. (f42"|) and write the nucleon propagator in the compact form: 



S(p) = S F (p) + Snip) 
S F (p)=\rpl + M*ip) 



Snip) = [r% + M*(p)J -^y5 (p° - E(p)) e (P F - |p|) , 



(45) 



(46) 
(47) 

(48) 



where E(p) is the single-particle energy, which is the solution of the transcendental equation: 



E(p) = \E*ip)-Z°ip) 



P0 =E(p) 



{p 2 [1 + £"(|p|, Eip))f + [M + S s (|p|, Eip))} 2 Y - £°(|p|, Eip)) . 



(49) 
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Note that we have assumed that the nucleon propagator has simple poles with unit residue. 
Within the approximation scheme we are working in this paper, this assumption is satisfied, 
as can be seen below. 

Following the scheme of neglecting the Feynman part of the nucleon propagator, Eq. 
we obtain for the energy density of nuclear matter the following expression: 



£w = y I d xT\ 

= 7 



■ w 



V.E.V. 



p p d 3 q q ■ q* + MM*(q) i cD 
u (2vr) 3 E*{q) 



(50) 



where V.E.V. means the vacuum expectation value of T$, and £^ and are the direct 
and exchange contributions, given by: 

rPf d 3 q M*(q)~ 



I-i 

2 



A. 

ml 



7 



cE 



Pf 



o (2tt) 3 E*(q) 
d 3 q 





f p p d 3 q 


V2 J ml 


1 Jo (2vr) 3 _ 



(51) 



2 1 Jo (2-KfE*(q) Jo (2n) 3 E*(k) 



d 3 k 



[glK{q-k) 



x 



1 

- - 1 

2 



[E{q)-E{k)fK{q-k) 



q*% + M*(q)M*(k) 



+2glK{q~k) 



1 

- - 1 

2 



[E{q)-E(k)] 2 \,{q-k) 



p t k* tl -2M*(q)M*(k)]} . 



(52) 



These expressions are very similar to the ones obtained in the Hartree-Fock approximation. 
Differences are contained in the fermion kinetic energy, the first term in Eq. (pPf), and in the 
factors f| — in Eqs. ( |51~D and (0). These differences arise because we are not using the 
nucleon field equation of motion. 

To implement the 6 expansion, we need to specify £ . We choose 



£o = ^(i7^-iW + ^ 



where 



M = M + n . 
The interpolated Walecka model is then given by 

C w (5) = C + 5 y-g u ih l y*il) + g*i>H + 
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(53) 



(54) 



(55) 



Notice that the 5 expansion technique could have also been applied to the meson fields 
explicitly. However, we have chosen to eliminate the meson fields in favor of the nucleon 
fields by means of Eqs. ( |36|) and (|37|) , therefore meson effects enter via the nucleon fields. 

Next we obtain the self-energy in perturbation theory, always neglecting vertex correc- 
tions and the Feynman part of the nucleon propagator. In zeroth order in 5, the nucleon 
self-energy, corresponding to the interpolated Lagrangian Eq. (53), is obviously zero: 



= . 

Therefore, the auxiliary quantities to be used in Eq. ( [46]) become: 

M*(p) = M = M + fi, 
P = P , 

(P°,P). 



(56) 



(57) 



The single-particle energy is simply given by: 



E(p)=E*(p) = E (p)= p 2 + M, 



(58) 



Using these in Eqs. Q5"ip and (0), we obtain for the zeroth order energy density the following 
expression: 

' p * d 3 q q 2 + MM 



£ (0) - -v / 
tw " 7 io 



(2tt)3 E (q) 
with the direct and exchange contributions given by: 



M0)D p (0)E 



(59) 



c(0)D 



mi 



7 



p F d 3 q M 
o (27r)3£ (g) 



1 ,\ //; 



7 



1 2 



o (2tt) 3 



(60) 



d 3 q 



d 3 k 



2 <y (2.).^) Jo w .^ (fc) {^(^)-W-(q-k) a ) 

i - l) - [£ (g) - J Eo(^)] 2 A (J ([E (g) - £o(A:)] 2 - (q - k) 2 )] [^(g)^(fc) 

i - l) - [EM E (k)] 2 



- q • k + M 2 j + 2^A w ([E (g) - £o(A:)] 2 - (q - k) 2 
x A^E^q) - -E'o(^)] 2 - (q - k) 2 )] [E (q)E (k) - q ■ k - 2M 2 ] } . 

Let us consider the direct term first. Application of the PMS to it: 



(61) 



15 



d/j dM dfx dM 
yields the following self-consistency condition for M : 



= , (62) 



d£$ d£% dM d£ w 



^ Pf d 3 g Mq 

M ~^:L wm- (63) 

This is exactly the same self-consistency condition for the effective nucleon mass obtained 
by means of the Hartree, or mean-field, approximation. 

Now, application of the PMS to the full energy density, which includes both direct 
and exchange contributions, leads to a nonlinear equation for /i, or equivalently for Mo, 
which is more complicated than the one of Eq. fl6"3"|). We do not present the expression 
here because it is rather lengthy and not very instructive. To investigate the size of the 
exchange corrections we carry out two sets of comparisons. In Figure 4 we compare the 
nucleon binding energy, obtained by using only the first and direct terms in Eq. (|59[) (solid 
line) and coupling constants fixed by fitting the binding energy and density of equilibrium 
nuclear matter, with the full binding energy, keeping the same coupling constants (dotted 
line). The value of the coupling constants are g 2 s = 91.64 and g% = 136.2. The masses 
used in all calculations are M = 939 MeV, m v = 783 MeV and m a = 550 MeV. We find 
that the exchange corrections coincide with those obtained in a relativistic Hartree-Fock 
calculation |TH|TBJ which we also show for comparison (long-dashed line). The dashed line 
shows the results obtained for the total binding energy (including both direct and exchange 
contributions) after renormalizing the model parameters to reproduce the bulk saturation 
properties of nuclear matter: g 2 s = 83.11 and g\ = 108.05. These coupling constants are 



the same used when renormalizing the relativistic Hartree-Fock calculation of Ref. [16 



Therefore, the PMS condition on the zeroth order energy density of the Walecka model is 
also equivalent to the usual Hartree-Fock solution. 

In Figure 5 we show the results for fjj cLS cL function of the Fermi momentum Pp ob- 
tained with the application of the PMS to the zeroth order energy density. The solid line 
corresponds to the first and direct terms only and the dashed one (almost unnoticeable) 
corresponds to the full energy density with the renormalized constants. In Figure 6 we 
compare the results for the effective nucleon mass in nuclear matter as a function of Pp 
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obtained from fi. In both figures, it is clear that the results with the exchange terms and 
renormalized constants coincide with the results obtained by using the direct terms only. 

We now consider the second-order contribution to the self-energy. The self-energy to 
second-order in delta is is given by: 



S f2, (p) 



d 4 q 



Tr [S (0) (<?)] 



g 2 J 2 r A 



mi 



+ ^ 2 j ^S^(q)A a (p-q)- ig 2 J 2 J ^^(^AUp-qh 



d 4 q 



7 M Tr[ 7 ^(°)(g) 



(27T) 



(64) 



where A CT and A w are given in Eq. (|3Bp, and again we have made use of the baryon cur- 
rent conservation. We evaluate this expression neglecting the Feynman part of the nucleon 
propagator, the term Sp(p) given by Eq. ((47]). Because of this, all integrals in Eq. (|64] ) are 
finite and can easily be evaluated; there is no need for renormalization. The first term in 
Eq. (|64D comes from the first order contribution in 5 and must be kept at second order. 



The results are very similar to the ones obtained with the Hartree-Fock approximation [16 



Since there are subtle differences, we write them explicitly below. Each component of the 
self-energy, E s , S°, and T, v , can be decomposed in a direct and an exchange part. The direct 
components are given by: 

g 2 S 2 Pf d 3 q Mo 



,s(2) 
J D 



,0(2) 
J D 



-6 (j, - 7 : 



mi 



(2vr) 3 £ (g) 



-1- 



g 2 J 2 r p F <P q 



mi 



(27T) 



E^ 2) = . 



(65) 

(66) 
(67) 



Notice that they are independent of energy and momentum. The exchange terms are given 
by: 



^f\p) 



,0(2) 



(P) 
(P) 



4n 2 p Jo 
1 

4:7T 2 p Jo 

1 

A-K 2 p 2 Jo 



dq q 



dq q 



Eo(q) 
1 



/L g 2 j 2 e a ( P ,q) + -g 2 j 2 e w (p,q) 



dq q 



q 



Eo(q) 



gj®c{p,q)+gj®u{p,q) 



where the functions Qi(p, q), &i(p, q),i = a,u, are defined by: 

A&q) + %pq 



Qi(p,q) = In 



Ai(p,q) - 2pq 



®i{p,q) = ^A(p, g)©<(p, ?) 



(68) 
(69) 
(70) 

(71) 
(72) 
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where 



Ai(p, q) = p 2 + q 2 + mj - \E{p) - E (q)f 



(73) 



The auxiliary quantities to be substituted into Eq. ( |46| ) are then given by: 



P =P 



E*{p) = p* 2 + M* 2 



(74) 



p*» = jf + E"(p) = p° + Sg 2) + Zf\p), p* 



and the single-particle energy is the solution of: 



E(p) = \E*(p) -£ 0(2) G») 



Po=E(p) 



p 2 [l + E v ^(\p\,E(p))f+[M + E s ^(\p\,E(p))] 2 y -E°^(\p\,E(p)). (75) 

We are in the position to calculate the energy density Initially we consider the direct 
terms only. The energy density is given by: 



~(2)D 



'W 



r p p d 3 q q • q* + MM* g 2 
7 Jo (2tt) 3 E*(q) + 2m u 



2 ,1 2 



3tt 



2 1 F 



2m 2 



7 



Pf d 3 q M* 
o (27r) 3 E*(g) 



(76) 



Application of the PMS to this yields again the familiar Hartree result, Eq. (p3|), with M 
and Eo(q) replaced respectively by M* and E*(q). The exchange terms are shown in Eq. (j52|) 
and we believe it is not necessary to rewrite them here. Also in this case, despite numerical 
imprecisions, the exchange corrections coincide with the usual Hartree-Fock solution, as can 
be seen in Figure 7. The behavior of M* as a function of the Fermi momentum at this 
order does not show any noticeable difference as compared with the zeroth order results. 
However, as can be seen in Figure 8 the behavior of /i as a function of Pp obtained with 
the application of the PMS to the full second order energy density (dashed line) is rather 
different from the one obtained when only the direct term is taken into account (solid line). 
It is very interesting to notice that this different behavior does not manifest itself either in 
the values of M* or of the binding energy. This is because the energy density £ is a very flat 
function of /i, as can be seen in Figure 9, where the energy density is shown for Pp = 1.19 
fm _1 . Recall that if one had an exact solution, the energy density would be independent of 



fi. The solid line is obtained without the inclusion of the exchange term (the PMS solution 
in this case is given by fi/M = —0.275) and the dashed line gives the full second order 
density energy (the PMS solution is [ijM = —0.35). This stability in the value of the energy 
density as a function of /i is very desirable and guarantees that even big changes in the value 
of fi will not affect physical quantities, as the binding energy for instance. 



IV. CONCLUSIONS AND FUTURE PERSPECTIVES 



In this paper we have utilized the optimized 5 expansion to study medium effects in 
two commonly used models in hadron and nuclear physics: the NJL model and the Walecka 
model. We have investigated an alternative way of fixing the arbitrary parameters introduced 
by the 5 expansion, by applying the PMS to the energy density of the system. 

The most important and concrete conclusion we can draw from this work is that when 
applying the PMS to the energy density of the NJL model we reproduce, at any order, the 
familiar Hartree-Fock solution for the dynamically generated mass, in the approximation of 
neglecting vertex corrections. In the case of Walecka model, we obtained results quantita- 
tively similar to the ones of the usual Hartree-Fock approximation, although the analytical 
expressions are not evidently equivalent. If one neglects the exchange term in the energy 
density then clearly the mean-field solution is reproduced at any order. It is also worth 
mentioning that, in the Walecka model, the energy density is a very flat function of \x and 
this guarantees that the PMS solution is indeed very stable. 

On the basis of our results, we believe that the optimized 5 expansion is a very robust 
nonperturbative approximation scheme. Compared with the Hartree-Fock approximation, 
the S expansion is very economical because of its perturbative nature. Once the reliability 
of the scheme has been established, one is ready to proceed to other interesting applications. 
These include vertex and, obviously, vacuum effects. The study of the vacuum in the Walecka 
model is an important issue since one needs to know the limits of applicability of such model 
to high densities and/or temperatures before quark and gluon degrees of freedom have to be 
invoked. Of course, one has to face renormalization problems when including the vacuum. 
Renormalization in a Hartree-Fock scheme is very complicated ]T7[ and one expects that 
this will be facilitated within the 5 expansion method. 
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FIGURES 



FIG. 1. Pp dependence of ji in the NJL model obtained with PMS on /„.. 

FIG. 2. Constituent quark mass (solid and dotted lines) and — < qq > : / 3 (dashed and 
dot-dashed lines) as a function of Pp. The solid and dashed lines are the PMS1 solution and 
the dotted and dot-dashed lines are the PMS2 solution. 

FIG. 3. Pp dependence of f w for the NJL model. The solid and dotted lines give respectively 
the PMS1 and PMS2 solutions with the same parameters. The dashed line gives the PMS2 solution 
with renormalized parameters. 

FIG. 4. Pp dependence of the binding energy of the Walecka model at zeroth order in 5. The 
solid line represents the first and direct terms of Eq. (59) only. The dotted and long-dashed lines 
give the full binding energy and the Hartree-Fock solution respectively, both determined with the 
same coupling constants used in the solid line solution. Finally the dashed line gives the full binding 
energy with the renormalized coupling constants. 

FIG. 5. Pp dependence of fi for the Walecka model at zeroth order in 5. The dashed line 
represents fi determined with the full self-energy, which is the sum of direct and exchange terms. 
The solid line represents \x determined without the exchange term. 

FIG. 6. Zeroth order nucleon effective mass M$ as a function of Pp. The solid curve is the 
result obtained without the exchange term and the dashed curve is the result using the full energy 
density. 

FIG. 7. Pp dependence of the binding energy of the Walecka model at second order in 5. The 
solid, dashed, dotted and long-dashed lines are the same as in fig. 4 

FIG. 8. Pp dependence of /i for the Walecka model at second order in 5. The solid and dashed 
lines are the same as in Figure 5. 
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FIG. 9. (i dependence of the energy density for the Walecka model at second order in 6, calcu- 
lated at Pp = 1.19 fm _1 . The solid line gives the solution when the exchange term is not included. 
The dashed line gives the full solution. 
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